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Abstract. 

The p and uj meson self-energy at finite temperature and baryon density have been analysed 
for an exhaustive set of mesonic and baryonic loops in the real time formulation of thermal field 
theory. The large enhancement of spectral strength below the nominal p mass is seen to cause 
a substantial enhancement in dilepton pair yield in this mass region. The integrated yield after 
space-time evolution using relativistic hydrodynamics with quark gluon plasma in the initial 
state leads to a very good agreement with the experimental data from In-In collisions obtained 
by the NA60 collaboration. 



1. Introduction 

Colliding heavy ions at ultra-relativistic energies is the only way to produce and study bulk 
properties of strongly interacting matter. Systematic efforts, both theoretical and experimental 
over the last few decades [H [2] have addressed various facets of the thermodynamics of the 
underlying theory — QCD. The spectra of hadrons emitted after freeze-out of the fireball 
produced in heavy ion collisions have provided us with a wealth of information. This includes 
the recent discovery made by studying the elliptic flow of hadrons that the quark gluon plasma 
(QGP) produced in Au-|-Au collisions at RHIC actually behaves as a strongly interacting fluid [3] 
as opposed to a asymptotically free gas of quarks and gluons. However, by virtue of the fact 
that electromagnetic probes (real photons and dileptons) are emitted all through the lifetime 
of the fireball coupled with their low rescattering probability make them penetrating probes 
capable of mapping the space-time history of the collision [3] . Lepton pairs with both invariant 
mass and transverse momentum information are in fact preferable to real photons. Large mass 
pairs produced by Drell-Yan process and from the decays of heavy quarkonia are emitted early 
whereas pairs with low invariant mass radiated from thermal hadronic matter and Dalitz decays 
of hadrons are produced late in the collision. The invariant mass spectra of dileptons thus carry 
time information as displayed explicitly in [5j using the invariant mass dependence of the elliptic 
flow of lepton pairs. 

The rate of production of thermal dileptons is proportional to the two-point correlator of 
vector currents ^ . In the low invariant mass region which is dominated by lepton pairs produced 
during the later stages of the collision, dilepton emission takes place due to the decay of vector 
mesons. Consequently, the spectral properties of vector mesons, the p meson in particular has 
been a subject of intense discussion [U HI [21 [S]. We find that it is only for the vr — vr loop that 



one calculates the thermal self-energy loop for the vector mesons directly. In the case of other 
loops typically involving one heavy and one light particle or both heavy particles one uses in 
general either the virial formula or the Lindhard function. 

The sources modifying the free propagation of a particle find a unified description in terms 
of contributions from the branch cuts of the self energy function. In addition to the unitary cut 
present already in vacuum, the thermal amplitude generates a new cut, the so called the Landau 
cut which provides the effect of collisions with the surrounding particles in the medium. This 
formalism was applied to obtain the p self-energy in hot mesonic [9j and baryonic [TU] matter 
considering an exhaustive set of one-loop diagrams. The framework of real time thermal field 
theory that we use, enables us to evaluate the imaginary part of the self-energy from the branch 
cuts for real and positive values of energy and momentum without having to resort to analytic 
continuation as in the imaginary time approach. To evaluate the baryonic loops we work with 
the full relativistic baryon propagator in which baryons and anti-baryons manifestly appear on 
an equal footing. Thus the contributions from all the singularities in the self-energy function 
including the distant ones coming from the unitary cut of the loops involving heavy baryons are 
also included. These are not considered in the Lindhard function approach but can contribute 
appreciably to the real part of the p meson self-energy as shown in the case of a A^A loop. 

The broadening of the vector meson spectral functions leads to an enhancement of lepton 
pair production in the invariant mass region below the p peak. The effect of the evolving matter 
is handled by relativistic hydrodynamics. The integrated yield is seen to agree very well with 
the NA60 data [I2j from In-In collisions at 17.3 AGeV. 

The article is organised as follows. We will begin with a short derivation of the the dilepton 
emission rate in terms of the current correlation function in section 2. The relation to the 
spectral function of vector mesons is specified in section 3. This will be followed by a discussion 
on p and oj self-energies in hot and dense matter in section 4. In section 5 a brief account of the 
space-time evolution and initial conditions will be provided followed by the dilepton invariant 
mass spectra. We will end with a summary in section 6. 



2. Dilepton emission rate and the current correlation function 

Let us consider an initial state \I) which goes to a final state \F) producing a lepton pair 
with momenta pi and p2 respectively. The dilepton multiplicity thermally averaged over initial 
states is given by 



^^'^ ' ' Z (2^)32^1(2^)32^2 ^ ' 

where Z = Tr[e~^^] and Cmt = e'ilJi{x)j^ipi{x)A'^{x) + eJ^(x)y4^(x) in which ipi{x) is the 
lepton field operator and J^{x) is the electromagnetic current of hadrons. Following 0113113] 
this expression can be put in the form 

_dN_ ^ e^^4f^^-PqoW> (n\(ni^a^ - n^^^ 



where, = J d'^x e'^'^'^{J^{x)J^{{))) p is the Fourier transform of the thermal expectation value 

of the two-point correlator of the hadronic currents and L{q ) = (1 H — ^)W1 We now 

define W^^ by interchanging the order of the currents getting W^^ = e~^'^° VF^ . Using these to 
define the commutator W^y = — = / d^x e*^'^([J^(x), Jy{^)])i3 we finally have 

'^^ - ^"'-L{q')fBEiqo)W^Aqo,q) (3) 
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for conserved hadronic currents. This expression appears in different forms in the hterature. 
Replacing fBE{qo)W^u in dSD by (i)W^J = / d-^x e^5-^( j/^(0))/3 yields the expression in [6], 
by (ii) 2ImW/;lfBE{qo) where W^^ = i ^ d'^x e*9-^0(i - t'){[J^{x), J!){0)])f^ yields the rate in [H] 
and by (iii) 21mWj^/(l + e^"?") where Wj^^ = ijd'^x e''?-^(TJ^(a;) j/^(0))/3 gives the rate in [E]. 

The rate given by eq. ([3]) is to leading order in electromagnetic interactions but exact to all 
orders in the strong coupling encoded in the current correlator W^i,. The in the denominator 
indicates the exchange of a single virtual photon and the Bose distribution implies the thermal 
weight of the source. 



3. Current correlator and the spectral function of vector mesons 

Thermal field theory is the appropriate framework to carry out perturbative calculations in 
the medium. In the real time version of this formalism two point functions assume a 2 x 2 
matrix structure on account of the shape of the contour in the complex time plane [16^ 117] . It 
is convenient to begin with the quantity, 

T^t = i J d'xe^'^-^T,j;:{x)4mp (4) 

where Tc denotes ordering along the contour and a, b are the thermal indices which take values 
1 and 2. This quantity can be diagonalised by means of a matrix U so that 




^J = crr ]u; c/= Kr , . „ = 3^ (5) 



where T^i, is an analytic function. Equating both sides it follows that this function is obtainable 
from any one (thermal) component of T^^. It is related e.g. to the 11-component as, 

ReTfj,^{qo, q) = Rer^^(gfo, q) ; lmr^,,(go, q) = coth(;^y^)Imr^^(go, q) ■ (6) 
Furthermore, T^j^ has the spectral representation [181 US] 

T (n a\- MgQ W^,{q'^,q) 

J 2tt q'^-qo- ir]e{qo) 



which immediately leads to 



Wf,„{qo,q) = 2e{qo)Im.T^^{qo,q) . (8) 



In the QGP where quarks and gluons are the relevant degrees of freedom, the time ordered 
correlation function T^l can be directly evaluated by writing the hadron current in terms of 
quarks of flavour / i.e. = X)/ ^/V'/Ta*^/- To leading order we obtain using relations and 

m, 

9''W,. = -'^j:e}il-'-^). (9) 

The rate in this case corresponds to dilepton production due to process qq ^ 'J* ^ To 
obtain the rate of dilepton production from hadronic interactions it is convenient to break up 
the quark current into parts with definite isospin 

Jjl = ^ {u-ff,u - d-f^d) + ^ {u-f^u + d-f^d) + ■■■ 

= J^ + Ju/^ + --- (10) 



where V and S denote iso- vector and iso-scalar currents and the dots denote currents comprising 
of quarks with strangeness and heavier flavours. These currents couple to individual hadrons 
as well as multiparticle states with the same quantum numbers and are usually labelled by the 
lightest meson in the corresponding channel |19j . We thus identify the isovector and isoscalar 
currents with the p and uj mesons respectively. Using eq. (fTO]l in eq. dH and neglecting possible 
mixing between the isospin states, we write 

T^l = Tlt' + Tli^^'''/^ + --- (11) 

where 

TjfJ''' = ^ / dV-^{TjP^{x)Jm)p (12) 
and similarly for the scalar current. The current commutator in the isospin basis follows as 

W^u = WP, + W^j9 + ... (13) 

The correlator of vector-isovector currents Wj^,^ have in fact been measured [20] in vacuum along 
with the axial-vector correlator by studying r decays into even and odd number of pions. The 
former is found to be dominated at lower energies by the prominent peak of the p meson followed 
by a continuum at high energies. The axial correlator, on the other hand, is characterised by 
the broad hump of the ai. The distinctly different shape in the two spectral densities is an 
experimental signature of the fact that chiral symmetry of QCD is dynamically broken by the 
ground state ^21j. It is expected that this symmetry may be restored at high temperature and/or 
density and will be signalled by a complete overlap of the vector and axial- vector correlators [2^ . 

In the medium, both the pole and the continuum structure of the correlation function gets 
modified [11[23]. We will first evaluate the modification of the pole part due to the self-energy 
of vector mesons in the following. Using Vector Meson Dominance the isovector and scalar 
currents are written in terms of dynamical field operators for the mesons allowing us to express 
the correlation function in terms of the exact (full) propagators of the vector mesons in the 
medium. Writing Jj^{x) = Fpmpp^{x) and J^{x) = 'iF^rriuiU: p{x) in eq. and using eq. ([8]) 
the current commutator becomes 

Wp, = 2e{qo)F^mllmD'^, + 2e{qo) F^mllmV;, + ■■■ (14) 

where D^i, is the diagonal element of the thermal propagator matrix which is a two point function 
of the fields of vector mesons and is diagonalisable as in ([5]). The exact propagator is obtained 
in terms of the in-medium self-energies using the Dyson equation. Following [24j this is given 
by 

2) ^(^q^ — -^M^ Qfiv/Q ^-j^g-j 

- 111?^ - Iit{q) q"^ - mj - q'^Ui{q) q'^mf, 

where P^i, and Q^^ are the transverse and longitudinal projections. The imaginary part is then 
put in eqs. (jl4p and then in eq. ([3|) to arrive at the dilepton emission rate [24] 

^4^4^ = -^^L{q^)fBE{qQ) FpmlAp{qo, q) + F^mlA^{qo, q) ^ (16) 

where e.g. Ap{= —g^^^lmD^ /?,) is given by 
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(g2-m2-^Ren«)2 + (^Imnf)2 (g2 _ ^2 - g2 ^ RgHf )2 + I^nf )2 

(17) 
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Figure 1. One- loop Feynman diagrams for p or uj self-energy involving mesons (first figure) 
and baryons (second and third figures). V stands for the p or w in the external line. In the 
internal lines, h stands for vr, uj, oi and hi mesons. For the baryonic loops, and B indicate 
respectively nucleon and baryonic internal lines. 




the sum running over all mesonic and baryonic loops. Thus, the dilepton emission rate in 
the present scenario actually boils down to the evaluation of the self energy graphs (shown in 
Fig- m). The self-energy is also a 2 x 2 matrix and is diagonalisable by the matrix U~^. The real 
and imaginary parts of the self energy function can then be obtained from the 11-component 
as [niEU] 

Ren^i/(^o,^) = Ren^;,(go,g) = Ren"^(go,g) 

Imn^;,(go, q) = e(9o)Imn^,.(go, q) = tanh(/3go/2)Imn^],(5o, q) (18) 

where denotes the retarded self-energy. 

As indicated earlier, coupling of the hadronic current to multiparticle states gives rise to a 
continuum structure in the current correlation function W^^ . Following Shuryak we take a 
parametrised form for this contribution and augment the dilepton emission rate with 

V=p,u 



where 



4Cont _ _ 1 _|_ n /'on^ 

P Stt V ^ TT y 1 + exp(a;o - go)/<5 ^ ' 

with ojQ = 1.3, 1.1 GeV for p, w and 5 = 0.2 for both p and oj. The continuum contribution for 
the OJ contains an additional factor of ^. 

4. Spectral function and self-energy of vector mesons 

Obtaining the in-medium spectral functions for the p and w mesons essentially involves the 
evaluation of the self-energies. For the p meson these have been evaluated recently for mesonic [9] 
and baryonic [IQj loops. In the following we summarise these results for the p followed by those 
for the cj meson. 

p meson 

The 11-component of the self-energy of vector mesons for loops containing mesons or baryons 
can be generically expressed as 

'^]}u{q) = i j j^L^,{k,q)A^\k,mk)A^\p,mp) . (21) 

where the term L^y{k,q) contains factors from the numerators of the two propagators in the 
loop as well as from the two vertices and 

A^^(A;, mk) = -tt, — + a^'Kibik^ - mhNi . (22) 
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Figure 2. Branch cuts of self-energy function in qo plane for fixed q given by vr/i loop. The 
quantities ^1,2,3 denote the end points of cuts discussed in the text : qi = \J {rrih + rriT^y + \q\'^, 
q2 = \/imh - m^y + and q3 = \q\- 

For bosons, a = 1 and Ni = n with n[ujk) = — ^ whereas for fermions, a = —1 and 

A^i = n-^9{k()) + n^9{—ko) with n±{uji:) = . The loop momentum p in (pT|) for the 

various cases are as shown in Fig. ([1]). 

We begin with mesonic loops. The p self-energy for four possible ir-h loops, where h = 
vr, w, hi, ai have been evaluated. The imaginary part of the retarded self-energy (fT8|l is given by. 



Iinn'^{qo,q} = j 



[L^i^iil + n^ + nh)5{qo - uj^r - ^^h) - {n^ - nh)5{qo - o;^ + ujh)} 

+Li^''{{n^ - nh)5{qQ + uj^, - uJh) - (1 + + nh)5{qQ + uj^ + w/i)}] . (23) 



where the Bose distribution functions n^^ = n{ujT^) with oj^^ = \J k'^ + and rih = n{oJh) with 



= y {q — ky + "T-h- = 1)2) are the values of L^^{kQ) for k^ = Uj^, —Uj^ respectively and 

the vertices used in L^^{kQ) have been obtained from the chiral Lagrangians which are specified 
in [9]. 

The regions, in which the four terms of eq. ()23p are non-vanishing, give rise to cuts in the 
self-energy function (Fig. [2]). These regions are controlled by the respective (5-functions. Thus, 
the first and the fourth terms are non-vanishing for q^ > {rrih + mT^)"^ , giving the unitary cut, 
while the second and the third are non- vanishing for q^ < (m/j — ttt-tt)^, giving the so-called 
Landau cut. The unitary cut arises from the states, which can communicate with the p. These 
states are, of course, the same as in vacuum, but, as we see above, the probabilities of their 
occurrence in the medium are modified by the distribution functions. On the other hand, the 
Landau cut appears only in medium and arises from scattering of p with particles present there. 
We note that this contribution appears as the first term in the virial expansion of the self-energy 
function. 

Integrating over the angle and restricting to the kinematic region qQ,q^ > the contribution 
to the imaginary part coming from the unitary cut is given by [S] 

1 f^^ 

Imn^'^ = y _ duj^Ll"" {I + n{uj^) + n{qQ - u^)} (24) 
and that from the Landau cut is given by 

+' 

Imll^'' = - / di^T^L^ {n{ujT,) - n(go + w^^)} (25) 



where a;^± = ^(go ± \^W^), a;^±' = §(-go T I^IW^) with = ^1-^^ and 




Figure 3. Left panel shows the imaginary part from baryonic loops and the right panel shows 
the total contribution from meson and baryon loops. 



The real part can be obtained from the imaginary part by a dispersion relation or can be 
evaluated directly from the graphs. 

Let us now turn to the baryonic loops in the p meson self-energy. Using p = k — q for the 
third diagram of Fig. [1] in eq. (|2ip the retarded self-energy is evaluated for NB loops including 
all spin one-half and three-half 4— star resonances listed by the Particle Data Group so that 
B stands for the A^*(1520), iV*(1650), iV*(1700), iV*(1720) A(1230), A*(1620), as well as the 
A^(940) itself. 

As before, the imaginary part can be evaluated from the discontinuities of the self-energy. 
However, the threshold for the unitary cut for the baryon loops being far away from the p pole 
we only consider the Landau part. The expression of the self-energy corresponding to the second 
and third diagrams of Fig. [1] can be obtained from one another by inverting the sign of q. Adding 
the contributions coming from the two diagrams, the imaginary part of the self energy is given 
by 

Imn^'" = "Yevrj^ ] + ^^A'[^r(-9){^+(^o+WAr)-n+(a;Ar)} + L^''(g){n_(go + t^Af)-?^-(wAr)}] 

(26) 

where = §(-go ± I^VF^v) with = ^1 - S% = q^ - m\ + m%. The factors 

(z = 1, 2) in eq. ([26]) are the values of L^'^iko) for = un, —ojn respectively and is obtained 
using gauge invariant interactions details of which are provided in [10] . 

The calculations described above treats the heavy mesons and baryon resonances in the 
narrow width approximation. We have included their widths by folding with their vacuum 
spectral functions as done in [25]. 

We now present the results of numerical evaluation. We plot in the left panel of Fig. [3] 
the imaginary part of the p self-energy for baryon loops as a function of the invariant mass 
= M for q = 300. The transverse (solid line) and longitudinal (dashed line) components 
Imllt and q^luiUi have been shown separately. The A^A^*(1520) loop makes the most significant 
contribution followed by the A^*(1720) and A(1700). On the right panel is plotted the spin- 
averaged p self-energy defined by 11 = ^(2nf +q'^Ili) showing contributions from the baryon and 
meson loops for two values of the baryonic chemical potential. The small positive contribution 
from the baryon loops to the real part is partly compensated by the negative contributions from 
the meson loops. The substantial baryon contribution at vanishing baryonic chemical potential 
reflects the importance of anti-baryons. 




Figure 4. The spectral function of the p meson for (left) different values of the temperature 
T and (right) different values of the baryonic chemical potential (//). 



We now turn to the spin averaged spectral function defined in eq. p7p . First, in the left panel 
of FigU] we plot the spectral function at fixed values of the baryonic chemical potential and 
three- momentum for various representative values of the temperature. We observe an increase 
of spectral strength at lower invariant masses resulting in broadening of the spectral function 
with increase in temperature. This is purely a Landau cut contribution from the baryonic loop 
arising from the scattering of the p from baryons in the medium. We then plot in the right panel 
of Fig. m the spectral function for various values of the baryonic chemical potential for a fixed 
temperature. For high values of ^ we observe an almost flattened spectral density of the p. 

4.2. uj meson 

The UJ self-energy is evaluated along similar lines [26) . The uj meson decays mostly into three 
pions. Assuming this decay to proceed via an intermediate p meson i.e. uj ^ pn ^ Stt, the 
dominant contribution to the uj self-energy in meson matter can be expressed as [27] 

(1) = ]^ dM' iKU M)]A,{M) (27) 

where Np = dM^Ap(M^), Ap being the vacuum spectral function of the p. Here 

^r'{itp) ('^' obtained by evaluating the first diagram of Fig. ([1]) with h = p. 

For UJ self-energy due to baryons we have evaluated NB loops where B = 
iV*(1440), A^*(1520),iV*(1535),iV*(1650), A^*(1720), A^(940). The calculation proceeds similarly 
as in case of the p. 

Left panel of Fig. ([5]) shows the transverse and longitudinal components of the imaginary part 
coming from the baryonic loops and the right panel shows the individual mesonic and baryonic 
loop contributions to the uj self-energy at two different chemical potentials. We see that the 
A'^*(1535) plays the dominating role mainly due to the strong coupling compared to the other 
baryonic resonances |27j. The spin averaged spectral function at different temperatures and 
chemical potentials are plotted in the left and right panels of Fig. ([6]) respectively. We observe 
a slight positive shift in the peak position. 

Having obtained the in-medium spectral functions of the two most important low-lying vector 
mesons, we are now in a position to evaluate the static rate of dilepton production using eq. (fT6]l . 
Integrating over the transverse momentum qt and rapidity y of the electron pairs we plot 
dR/dM'^ vs M in Fig. [7] for T=175 MeV. Because of the kinematical factors multiplying the p 





Figure 6. The spectral function of the w for different values of T (left panel) and (right 
panel) . 



spectral function the broadening appears magnified in the dilepton emission rate. A significant 
enhancement is seen in the low mass lepton production rate due to baryonic loops over and above 
the mesonic ones shown by the dot-dashed line. The substantial contribution from baryonic 
loops even for vanishing chemical potential points to the important role played by antibaryons 
in thermal equilibrium in systems created at RHIC and LHC energies. 

5. Space time evolution and dilepton spectra in In-In collisions 

Since dileptons are produced at all stages of the collision it is necessary to integrate the emission 
rates over the space-time volume from creation to freeze-out. We assume that quark gluon 
plasma having a temperature Tj is produced at an initial time Tj. Hydrodynamic expansion and 
cooling follows up to a temperature Tc where QGP undergoes a transition to hadronic matter. 
Subsequent cooling leads to freeze-out of the fluid element into observable hadrons. In the 
present work the fireball is taken to undergo an azimuthally symmetric transverse expansion 
along with a boost invariant longitudinal expansion p8]. The local temperature of the fluid 
element and the associated flow velocity as a function of the radial coordinate and proper 
time is obtained by solving the the energy momentum conservation equation 5^ T^^ = where 
T^^ = {e + P)u^vy + g'^'^P is the energy momentum tensor for ideal fluid. This set of equations 
are closed with the Equation of State (EoS); typically a functional relation between the pressure 
P and the energy density e. It is a crucial input which essentially controls the profile of expansion 




Figure 7. The lepton pair emission rate at T = 175 MeV with and without baryon (B) loops 
in addition to the meson (M) loops. 



of the fireball. 

The initial temperature is constrained by the experimentally measured hadron multiplicity 
through entropy conservation [29] , 

27r^ dN 
45C(3)7ri?^4afe^"d^ 



Thb^)ri = o2.. i—ib^)) (28) 



f{T)sg + (1 - fiT))sh with /(T) = -(1 + tanh(— ^)) . (29) 



where {dN / dy{hm)) is the hadron (predominantly pions) multiplicity for a given centrality class 
with maximum impact parameter 6^, R± is the transverse dimension of the system and is 
the degeneracy of the system created. The initial radial velocity, Vr{Ti,r) and energy density, 

e{Ti,r) profiles are taken as [30], Vriji^r) = and e{Ti,r) = eo/(e « +1) where the surface 
thickness, 5 = 0.5 fm. In the present work we assume Tc = 175 MeV [3T]. In a quark gluon 
plasma to hadronic matter transition scenario, we use the bag model EoS for the QGP phase and 
all resonances with mass < 2.5 GeV for the hadronic gas. The transition region is parametrized 
as [32] 

1 T — T 

-(l+tanh(^ " 

where Sg (sh) is the entropy density of the quark (hadronic) phase at T^. The value of the 
parameter T can be varied to make the transition strong first order or continuous. We take 
r = 20 MeV in this work. 

The ratios of various hadrons measured experimentally at different ^snn indicate that the 
system formed in heavy ion collisions chemically decouple at a temperature (T^h) which is 
higher than the temperature for kinetic freeze-out (Tj-) determined by the transverse spectra 
of hadrons [33]. Therefore, the system remains out of chemical equilibrium from T^h to Tj. The 
chemical non-equilibration affects the dilepton yields through (a) the emission rate through the 
phase space factor and (b) the space-time evolution of the matter via the equation of state. 
The value of the chemical potential and its inclusion in the EoS has been taken into account 
following p3] . 

Finally, we have obtained the dimuon yield {dN/dM) in In-In collisions at SPS at a center of 
mass energy of 17.3 AGeV. The initial energy density is taken as 4.5 GeV/fm^ corresponding to 
a thermalisation time = 0.7 fm. We take the QGP to hadronic matter transition temperature 
Tc = 175 MeV and the freeze-out temperature Ty =120 MeV which can reproduce the slope 
of the hadronic spectra measured by the NA60 Collaboration. In Fig. ([8]) we have shown the 
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Figure 8. Dilepton invariant mass spectra for different p^-bins compared with the NA60 data. 



invariant mass spectra for different transverse momentum (px) windows. The theoretical curves 
agree quite weU with the experimental data [12] for all the pT ranges. The strong enhancement 
in the low M domain is clearly due to the large broadening of the p (and w) in the thermal 
medium which comes entirely from the Landau cut in the self-energy diagrams. In the last panel 
we also plot for comparison the spectra calculated in [33] where the self-energy due to baryons 
has been evaluated following the approach of [M]- It is seen that this approach depicted by 
the dashed curve does not produce the required enhancement to explain the data in the range 
0.35 <M< 0.65 GeV. 



6. Summary 

The self-energy of p and uj mesons have been computed in nuclear matter at finite temperature 
and baryon density. Loop graphs involving mesons, nucleons and 4-star A^* and A resonances 
up to spin 3/2 were calculated using gauge invariant interactions in the framework of real 
time thermal field theory to obtain the correct relativistic expressions for the self-energy. The 
singularities in the complex energy plane were analysed and the imaginary part obtained from the 
Landau cut contribution. Results for the real and imaginary parts at non-zero three-momenta 
for various values of temperature and baryonic chemical potential were shown for the individual 
loop graphs. The spectral function of the p was observed to undergo a significant modification 
at and below the nominal rho mass which was seen to bring about a large enhancement of 
lepton pair yield in this region. After a space-time evolution using relativistic hydrodynamics, 
the invariant mass spectra for various pT windows was found to be in very good agreement with 
the experimental data obtained in In-In collisions at 17.3 AGeV. 

References 

[I] "The Physics Of The Quark-Gluon Plasma," eds. S. Sarkar, H. Satz and B. Sinha, Lect. Notes Phys. 785, 1 

(2010). 

[2] "The CBM Physics Book," eds. B. Friman et al, Lect. Notes Phys. 814, 1 (2011). 
[3] L. P. Csernai, J. I. Kapusta and L. D. McLerran, Phys. Rev. Lett. 97 (2006) 152303. 
[4] J. Alam, S. Sarkar, P. Roy, T. Hatsuda, B. Sinha, Ann. Phys. 286, 159 (2000). 

[5] P. Mohanty, V. Roy, S. Ghosh, S. K. Das, B. Mohanty, S. Sarkar, J. Alam and A. Chaudhuri, Phys. Rev. C 

85 (2012) 031903. 
[6] L. D. McLerran and T. Toimela, Phys. Rev. D 31, 545 (1985). 
[7] R. Rapp and J. Wambach, Adv. Nucl. Phys. 25, 1 (2000) 
[8] S. Leupold, V. Metag and U. Mosel, Int. J. Mod. Phys. E 19, 147 (2010) 
[9] S. Ghosh, S. Sarkar and S. Mallik, Eur. Phys. J. C 70, 251 (2010) 
[10] S. Ghosh and S. Sarkar, Nucl. Phys. A 870-871, 94 (2011) 

[II] S. Ghosh, S. Sarkar and S. Mallik, Phys. Rev. C 83, 018201 (2011). 

[12] R. Arnaldi et al for the NA60 collaboration, Phys. Rev. Lett. 100 022302 (2008); Eur. Phys. J. C 61, 711 
(2009); S. Damjanovic et al for the NA60 Collaboration, J. Phys. G: Nucl. Part. Phys. 35, 104036 (2008). 
[13] H. A. Weldon, Phys. Rev. D 42, 2384 (1990). 

[14] R. Rapp, J. Wambach and H. van Hees, arXiv:0901.3289 [hep-ph]. 
[15] J. V. Steele, H. Yamagishi and I. Zahed, Phys. Lett. B 384 (1996) 255. 
[16] S. Mallik and S. Sarkar, Eur. Phys. J. C 61, 489 (2009) 

[17] M. Le Bellac, Thermal Field Theory (Cambridge University Press, Cambridge, 1996). 
[18] S. Sarkar, B. K. Patra, V. J. Menon and S. Mallik, Indian J. Phys. 76A, 385 (2002). 
[19] E.V. Shuryak, Rev. Mod. Phys. 65, 1 (1993) 

[20] R. Barate et al. [ALEPH CoUaboration], Eur. Phys. J. C 4, 409 (1998), K. Ackerstaff et al. [OPAL 

Collaboration], Eur. Phys. J. C 7, 571 (1999) 
[21] S. Sarkar, Nucl. Phys. A 862-863 (2011) 13. 
[22] J. I. Kapusta and E. V. Shuryak, Phys. Rev. D 49, 4694 (1994). 
[23] S. Mallik and S. Sarkar, Eur. Phys. J. C 25, 445 (2002). 
[24] S. Ghosh, S. Sarkar, J. Alam Eur. Phys. J. C 71, 9 (2011). 
[25] S. Sarkar, E. Oset and M. J. Vicente Vacas, Nucl. Phys. A 750, 294 (2005) 
[26] S. Ghosh, S. Mitra and S. Sarkar, (under preparation) 

[27] P. Muehlich, V. Shklyar, S. Leupold, U. Mosel, M. Post, Nucl. Phys. A 780, 187 (2006) 

[28] J. D. Bjorken, Phys. Rev. D 27, 140 (1983). 

[29] R. C. Hwa and K. Kajantie, Phys. Rev. D 32,1109 (1985). 

[30] H. von Gersdorff, M. Kataja, L. D. McLerran and P. V. Ruskanen, Phys. Rev. D 34 794(1986). 
[31] S. Borsanyi et al, Jour. High Eener. Phys. 1011, 077 (2010). 
[32] M. Asakawa and T. Hatsuda Phys. Rev. D 55 4488 (1997). 

[33] P. Braun-Munzinger, J. Stachel, J. P. Wessels and N. Xu, Phys. Lett. B 344 43 (1995); 365 1 (1996). 
[34] T. Hirano and K. Tsuda, Phys. Rev. C 66 054905 (2002). 

[35] J. K. Nayak, J. e. Alam, T. Hirano, S. Sarkar and B. Sinha, arXiv:0902.0446 [nucl-th]. 
[36] V. L. Eletsky, M. Belkacem, P. J. Ellis and J. I. Kapusta, Phys. Rev. C 64 (2001) 035202 



reserved for figure 



